Dynamics of gravitational field within a wave front and thermodynamics of black holes 
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I. INTRODUCTION 

Evolution of gravitational field within a finite tube 
with a timelike boundary was derived using a consistent, 
Hamiltonian formulation in [1] and then reformulated in 
[2] . Here we extend this description to the case of a null- 
like boundary or a wave front, i.e., a three-dimensional 
submanifold whose internal metric is degenerate. Re- 
stricting our result to the special case of wave fronts, 
namely to non-expanding horizons, we obtain a general- 
ization of the first law of thermodynamics for black holes 
as a simple consequence. 

Contrary to the Iyer and Wald approach (see [3]), 
no assumption about stationarity (existence of a Killing 
field) is necessary here. Such an assumption finally re- 
duces our formula to the standard first law. 

In many presentations of the Hamiltonian field theory 
(cf. [4]) boundary problems are neglected. Consequently, 
all the surface integrals arising from the integration by 
parts are assumed to vanish. Here we use the formulation 
proposed in [5] and [6] , where the field boundary data are 
treated on the same footing as the Cauchy data. This is 
the only way to obtain a mathematically consistent (in- 
finite dimensional) Hamiltonian description of any field 
theory if the boundary of the space volume taken into 
account is non-trivial. 

To illustrate our approach, consider as an example the 
linear theory of an elastic, finite string. Field configura- 
tion of the string is described by its displacement func- 
tion: K x [a,b] 3 (t,x) — > <p(t,x) e R, fulfilling the wave 
equation: 
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c 2 dt 2 ^ dx 2 ^ 



(1.1) 



Here, velocity "c" is a combination of the string's proper 
density (per unit length) and its elasticity coefficient. 
Passing to appropriate time and length units, we may 
always put c — 1. Equation (1.1) may be derived from 



the Lagrangian density 
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L = --v^det 



«T(t^)(M 



(1.2) 



where /i, v — 0,1 and (a; , a; 1 ) = (i, x), g M „ = 
diag(— 1, +1), "dot" denotes the time derivative and 
"prime" denotes the space derivative. A convenient way 
to encode the entire information about the dynamics of 
the string is to write it in a form of the following equation: 

SL(<p, d vV ) = d M (p^cp) = (V) 6<p + P »5 (o» , (1.3) 

equivalent to the Euler-Lagrange equations. Indeed, con- 
dition (1.3) is satisfied if and only if the volume part of 
the variation of L (normally present on the right-hand 
side) vanishes. The above equation (which we refer to, 
as the Lagrangian generating formula of the dynamics) 
has a beautiful symplectic interpretation (see [5] or [6]) 
as a definition of a Lagrangian submanifold of physically 
admissible states within the symplectic space of first jets 
of sections of the state bundle. Here, (ip,d v (p,d IJ ,p tl ,p 1 ') 
are local canonical coordinates in this symplectic space. 
Without going into deep, mathematical details, formula 
(1.3) may be simply read as the following condition im- 
posed on the string configuration: the "response param- 
eters" (d^p^jp^) must be equal to partial derivatives of 
the Lagrangian L with respect to the corresponding "con- 
trol parameters" (ip,d v tp). Hence, we have the following 
dynamical equations of the theory: 

1 . definition of the canonical momenta: 



• kinetic momentum 



P 



dL 



d(d ip) 
• stress density 



dof = <f>, 



dL 



P 



d{jh<p) 



= -d-^cp = -Lp 1 , 
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2. Euler-Lagrange equation, equivalent to (1.1): 



For purposes of the Hamiltonian description of the theory 
we introduce the following notation: 



7r := p° ; it 1 - := p 1 . 



(1.4) 



Integrating infinitesimal generating formula (1.3) over 
the entire string [a, b] we obtain the spatially finite gener- 
ating formula (it is still infinitesimal with respect to the 
time variable): 

5 [ L= [ (-kSp + nSip) + [^Sip]^ . (1.5) 

J a J a 

Hamiltonian description of the same dynamics is ob- 
tained via Lcgcndrc transformation between tt and tp, 
putting irSip — 5(irip) — ipSn: 

r b h 

-SH= (ir5ip - ipSn) + [ir ± 5<p] a , (1.6) 

J a 



with 



n 



■= f\n<p-L)= 1 - ftf + tff). (1.7) 

J a ^ J a 



Equation (1.6) acquires an infinitely-dimensional, Hamil- 
tonian meaning: 



7T = — 



5H 

Sip ' 



. _ 5H 

^ Sir ' 



(1.8) 



as soon as the boundary term in (1.6) is killed. This 
may be done in many ways, by imposing appropriate 
boundary conditions. Physically, this corresponds to 
a choice of a specific device controlling the behaviour 
of the extremal points of the string. Mathematically, 
such a choice implies a specific self adjoint extension 
of the second order differential operator on the right- 
hand side of the dynamical equation (1.1). This makes 
the dynamics uniquely defined. As an example con- 
sider two such choices: Dirichlet conditions and the Neu- 
mann conditions. In the Dirichlet mode we restrict our- 
selves to an infinitely dimensional phase space of ini- 
tial data (<p,tt), defined on [a, b] and fulfilling conditions: 
ip(a) = A, ip(b) = B. Within this phase space we have 
Sip (a) = Sip (a) = and equations (1.8) hold. 

Instead of controlling the string configuration tp at 
the boundary, we may control its stress by applying 
an appropriate force F. This leads to the Neumann 
control mode: tt (a) = F\ c f t , 7r (6) = fright- Conse- 
quently, Sir 1 - vanishes at the boundary. Performing Leg- 
endre transformation between ip and tt 1 - at the bound- 
ary: n Sip = S^^ip) — (pS'K ± , we obtain again a legitime 
Hamiltonian system, defined in a phase space which is 



completely different from the previous one and with a 
new Hamiltonian 7i playing role of free energy: 



—SH = I (ir Sip — ip Stt) — [<p 

J a 



Stt j 



(1.9) 



where 



H:=H+[ip ^] h a = H-[ip ip't 

=\ fy-wf-ipp"). (i.io) 

Again, boundary term in (1.9) vanishes due to Neumann 
conditions and the field dynamics reduces to (1.8). 

Consider now the subspace of static solutions: re — 
= p. Due to (1.8), the functional derivative of the 
H vanishes at those points and, whence, Hamiltonian 
formula (1.6) describes only the virtual work performed 
by the configuration controlling device at the boundary: 



sn = -[^sip] a . 



(i.ii) 



But, due to (1.7), H is manifestly convex. This implies 
that every static solution corresponds to the minimal 
value of the Hamiltonian in the phase space defined by 
the Dirichlet conditions. Due to equation (1.1) and to 
boundary conditions, such a solution is given by: tt = 
and ip(x) = A + (x — a) B b Z^ ■ Inserting this value into 
(1.7) we obtain the following "Penrose- like inequality": 



{B-Af 



<H, 



(1.12) 



analogous to the gravitational Penrose inequality relating 
the energy carried by Cauchy data outside of a horizon S 
and the energy of a black hole corresponding to the same 
value of appropriate boundary data on S. 

In the Neumann mode, the new Hamiltonian H is ob- 
viously non-convex. It is easy to check, that the free en- 
ergy (1.10) is unbounded neither from below nor from 
above and possesses no stationary points as soon as 
F\eft — bright ^ 0. There is, therefore, no "Penrose-like" 
inequality in this mode. 

In [1] the dynamics of the gravitational field within a 
timelike world tube S was analyzed in a similar way. For 
this purpose the so called "affine variational principle" 
was used, where the Lagrangian function depends on the 
Ricci tensor of a spacetime connection V. In this picture, 
the metric tensor g arises only in the Hamiltonian for- 
mulation as the momentum canonically conjugate to T. 
Later, it was proved in [2] that the Hamiltonian dynamics 
obtained this way is universal and does not depend upon 
a specific variational formulation we start with (actually, 
it can be derived from field equations only, without any 
use of variational principles, the existence of them being 
a consequence of the "reciprocity" of Einstein equations 
- see [5] and [6]). On the contrary, the Hamiltonian pic- 
ture is very sensitive to the method of controlling the 



3 



boundary data. A list of natural control modes, lead- 
ing to different "quasilocal Hamiltonians" , is given in [2] . 
Each of them is related with a specific choice of control 
variables at the boundary. The "true mass" , which tends 
to the ADM mass when shifting the boundary to infin- 
ity, is one of them (see also an analysis of the linearized 
theory [7]). 

The aim of the present paper is to generalize above 
results to the case when the boundary S is a wave front 
(a three-dimensional submanifold of spacetime M whose 
internal three-metric g a b is degenerate). This way we ob- 
tain a general Hamiltonian formula for the gravitational 
field dynamics within (or outside) a wave front, which is 
very much analogous to formula (1.6) for the string the- 
ory. As a byproduct, assuming that the wave front S is 
very special, namely is a non-expanding horizon, we ob- 
tain a generalization of the first law of thermodynamics 
for black holes (see formula (4.1)). Some of results have 
been already published in [8] . 



II. DYNAMICS OF THE GRAVITATIONAL 
FIELD WITHIN A NULL HYPERSURFACE 

Consider gravitational field dynamics inside a null hy- 
persurface S: 

s>0 s<0 s<0 s>0 




Parameter s = ±1 labels two possible situations: an 
expanding or a shrinking wave front (if S is a horizon, 
these correspond to a black hole or a white hole case). 
To simplify notation we use coordinates n = 0, 1, 2, 3, 
adapted to the above situation: x° = t is constant on a 
chosen family of (spacelike) Cauchy surfaces S t whereas 
x 3 is constant on the boundary S. This does not mean 
that x 3 is null-like everywhere, but only on S. We 
may imagine that x 3 = r coincides, far away from S, 
with the spacelike radial coordinate. Consider the three- 
dimensional volume V <zY, ta defined as the set of those 
points of S to which are situated inside S. Coordinates 
x A , A — 1,2, are "angular" coordinates on the two- 
surface dV = V C\ S whose topology is assumed to be 
that of a two-sphere. Finally, x k , k = 1,2,3, are spatial 
coordinates on the Cauchy surfaces {x° = const.} and 
x a , a = 0, 1, 2, are coordinates on 5*. We stress, however, 
that our results are coordinate independent and will be 
expressed in terms of relations between geometric objects 
defined on V and S. 

In Appendix B we prove the following identity ful- 
filled by any one-parameter family of solutions of Einstein 
equations ("variation" operator 5 may be understood as 
a derivative with respect to this parameter and "dot" 



denotes the time derivative): 
-8H J (P kl 5g kl - g kl 8P kl ) 

87r JdV 

+ 77T / {X AB &9ab - 2 (w SA° - A A Sw A )) , (2.1) 

J-D 71 " JdV 

where 

K=±- [ G"o + £ I U=-£-f A/, (2.2) 

87T Jy »7T J d y »7T J gy 

and P kl denotes external curvature of the Cauchy sur- 
face, written in the ADM form (cf. [9]). Moreover, 
A = \J det gAB is the two-dimensional volume form on dV 
and a — — \ log |g 00 | = log N, where N is the lapse func- 
tion. To define the remaining objects we must choose 
a null (time oriented) field K tangent to S. It is not 
unique, since fK is also a null (time oriented) field for 
any (positive) function / on S. For purposes of the 
Hamiltonian formula (2.1) we always choose the normal- 
ization compatible with the (3+l)-dccomposition used 
here: < K, dx° > = 1. Hence, K = d - n A d A . On 
the contrary, the vector-density A a = XK a = (A, —An' 4 ) 
is uniquely defined on S and does not depend upon a 
particular choice of K. Now, we define 1 

Lb ■= - g(db, Vaif) = -^£ K g a b, (2.3) 

w a :=- <V a K,dx° >, (2.4) 

where g a b is the induced (degenerate) metric on S. Be- 
cause of the identity: l a bK a = 0, the null mean curva- 
ture: I — g AB lAB may be defined (it is often denoted by 
9 - see [12], [11]), where by g AB we denote the inverse 
two-metric. 

The volume term in (2.2) vanishes due to constraint 
equations 2 g° v = 0. G° is often denoted by NH + N k H k 
(see e.g., [13]), where H is the scalar ("Hamiltonian") 
constraint and H k are the vector ("momentum") con- 
straints, TV and TV* are the lapse and the shift functions. 
Constraint equations H = and H k = imply vanishing 

In the Appendices B and C we give two independent 
proofs of the identity (2.1). The first one is analogous to 
the transition from formula (1.3) to formula (1.6). For 
this purpose we use Einstein equations written analo- 
gously to (1.3) (cf. [1]): 

5L = 8 K (w^8A%) , (2.5) 



1 This is the (2+l)-dccomposition of the extrinsic curvature 
Q a b(K) defined in [10] and [11]. 

2 In the presence of matter the volume term equals Q° — 8ttT°q 
and also vanishes due to constraint equations. 
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and A\ v 



flis 



°(p l v)k- 



where tt"" := ^Vlffl <T : 

Integrating (2.5) over a volume V and using metric con- 
straints for the connection T, we directly prove (2.1). 

However, in Appendix C, an indirect proof is also pro- 
vided, based on a limiting procedure, when a family S e 
of timelike surfaces tends to a lightlikc surface S. It is 
shown that the non-degenerate formula derived in [1] and 
[2] gives (2.1) as a limiting case for e — > 0. 

The last term in (2.1) may be written in the following 
way: 

-A A 5w A = Xn A SwA = n A SWA - u a waS\ 



III. DYNAMICS OF GRAVITATIONAL FIELD 
OUTSIDE OF A NULL SURFACE 

Consider now dynamics of the gravitational field out- 
side of a wave front S~ . We first add an external, time- 
like (non-degenerate) boundary S + and the situation is 
illustrated by the following figure: 



s<0 



s>0 



where Wa '■= Xwa and n A := g AB goB- Denoting n := 
ti a wa — wo = —K a w a we finally obtain the following 
generating formula: 




s>0 



s<0 



dV+ 




dV+ 



SH=^ j (P kl 6g kl -9kiSP kl ) 
8?r JdV 

+ lkj dv ( XlAB6 9 AB + 2 ( kSX nA6W *)) ■ -SH = -SH + SH 



(2.6a) 
(2.6b) 



where dV+ = V n S+, and dV~ = V n S~ . Be- 
cause dV~ enters with negative orientation, we have: 
J gv = j gv+ — jgy- ■ Integrating again Einstein equa- 
tions written in the form (2.5), over V, we use techniques 
derived in [1] and [2] to handle surface integrals over the 
timclike surface S + . To handle surface integrals over S~ 
we use our formula (2.1). This way we obtain: 



(2.6c) 

It is easy to prove that the integral lines of K are null 
geodesies. This means that K a V a K is always propor- 
tional to K. Hence, quantity k (traditionally called a 
"surface gravity" on S) fulfills equation K a \7 a K = kK, 
which may be used as its alternative definition. We stress 
that its value does not correspond to any intrinsic prop- 
erty of the surface S, but depends upon a choice of the 
null field K on S (i.e., upon a (3+l)-dccomposition of 
spacetime). However, in a special case of a black hole 
thermodynamics, there is a privileged choice of K, com- 
patible with the Killing field of the stationary solution 
and its normalization to unity at infinity. In this case k 
is an intrinsic property of the hole and the above formula 
provides, as will be seen later, the so called first law of 
black hole thermodynamics. 

We stress that the symplectic structure of gravitational 
Cauchy data is given here by the two first terms on the 
right-hand side of (2.6). Neglecting the second (surface) 
integral, the symplectic form would not be gauge invari- 
ant with respect to spacetime diffeomorphisms (see [2]). 
The sum of these two terms plays, therefore, role of the 
integral over the string interval [a, b] in formula (1.6). 
Most authors analyzing these problems take only the first 
(volume) integral as the symplectic form, which makes 
the entire approach gauge-dependent. 

The last integral (2.6c) is responsible for the control 
of five components of the boundary data: the two-metric 
qab on dV and the "curvature" Wa- Assuming the null- 
like character of S we already control the sixth parame- 
ter: g 33 = (see formula (A2)). This corresponds to the 
general observation (cf. [2]) that we must always control 
four gauge parameters of the boundary S, together with 
boundary data of the two "true degrees of freedom" of 
the gravitational field. 



lLJ v ( pklS ^-^ Spkl ) 
If (x 5a - a SX) + ^[ (XSa-aSX) 

>7T JdV+ V / 87T Jgy- 



77T- / Q ab S 9ab 

lOTT J d y+ 



/ {XI SgAB 



2(w SA° - A A Sw A )) , (3.1) 



where a is the "hyperbolic angle" between V and S + , 
whereas Q ab is the external curvature of S + written in 
the ADM form (cf. [2]). The contribution H + to the 
total Hamiltonian from the external boundary is written 
here in the form of a "free energy" proposed in [2] : 



i-\ Q°o-E Q , 



(3.2) 



1 

8 71 " Jdv+ 

where the additive gauge Eg is chosen in such a way that 
the entire quantity vanishes if dV + is a round sphere in 
a flat space. The internal contribution TC~ to the energy 
is given by formula (2.2) with dV replaced by dV~ . It 
was proved in [2] that shifting the external boundary to 
space infinity: dV + — > oo, the external energy H + gives 
the ADM mass, which we denote by M, whereas the 
remaining surface integrals over dV + vanish. Using this 
procedure we obtain the following generating formula for 
the field dynamics outside of an arbitrary wave front S~ 
in an asymptotically flat spacetime: 

-6M - SH- = ^ J (P kl Sg kl - g kl 8P kl ) 

+ — / (XSa - aSX) 
8lT Jov- 



+ 



-f- / (Xl AB SgAB + 2( K 5X-n A SW A )). (3.3) 

lOTT JdV- 
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IV. BLACK HOLE THERMODYNAMICS 

In this Section we apply the above result to the situa- 
tion, when the wave front S~ is a non-expanding horizon, 
i.e., I = (see [11]). In this case the "internal energy" 
?i~ given by formula (2.2) vanishes. Moreover, Einstein 
equations imply l AB = and definition (2.4) of w a re- 
duces to: V a K = —w a K (see [14]). Hence, we obtain the 
following generating formula for the black hole dynamics: 

SM = ^j v (P kl S 9 ki - 9kiSP M ) 

+ F~ / i^ Sa ~ hSX ) 
8lT Jov- 

+ ± { ( K S\-n A 6W A ), (4.1) 
71 " Jdv- 

where s = 1 for a white hole, and s = — 1 for a black 
hole. 

The so called "black hole thermodynamics" consists 
in restricting the above analysis only to stationary situ- 
ations. By stationarity we understand the existence of 
a timelike symmetry (Killing) vector field outside of the 
horizon. If such a field exists, we may always choose 
a coordinate system such that the Killing field becomes 
oItt and all the time derivatives (dots) vanish. Hence, 
formula (4.1) reduces to: 

SM = f (nSX - n A 5W A ) ■ (4.2) 

87r JdV- 

Moreover, we assume that is tangent to S. If this 
was not the case, we would have had a one-parameter 
family of horizons. Such phenomenon corresponds to the 
Kundt's class of metrics (see e.g., [15]). The known met- 
rics of this class are not asymptotically flat. However, we 
do not know whether or not this is a universal property 
and we exclude such a pathology by the above assump- 
tion. 

We have shown in [14] that there is a canonical affinc 
fibration n : S — > B over a base manifold B, whose topol- 
ogy is assumed to be that of a sphere S 2 . The affinc 
structure of the fibers is implied by the fact that they are 
null geodesic lines in M. Identity: — 2l a b = £-Kgab = 
implies that the metric g on S may be projected onto 
the base manifold B, which acquires a Rlemannian two- 
metric tensor Kab- The degenerate metric g a b on the 
manifold S is simply the pull back of Hab from B to S: 
g = n*h. 

The quantity w a is not an intrinsic property of the 
surface itself, but depends upon a choice of the null field 
K on S. Indeed, if K = cxp(—j)K then w a — w a + d a j. 
In particular, there are on S vector fields K such that 
K a \7 a K = and, consequently, k = 0. They correspond 
to the affine parametrization of the fibers of tt : S — > B. 

In case of a black hole, there is a privileged field K, 
compatible with the timelike symmetry of the solution, 
which is normalized to unity at infinity. This way the 



quantities k and wa in formula (4.2) become uniquely 
defined. 

We have, therefore, two symmetry fields of the metric 
g a b on S: 8q and K. Due to normalization chosen above, 
we have < do — K,dx° >= 0. Hence, the field n := 
do — K = n A dA is purely spacelike and projects on B. 
Moreover, it is a symmetry field of the Riemannian two- 
metric Hab- 

Because the conformal structure of Iiab is always iso- 
morphic to the conformal structure of the unit sphere 
5 2 , we are free to choose a coordinate system in which 
fiAB = fh-AB (and Hab denotes the standard unit two- 
sphere metric). The field n is, therefore, the symmetry 
field of this conformal structure. Consequently, n be- 
longs to the six-dimensional space of conformal fields on 
the two-sphere. Using remaining gauge freedom, we may 
choose angular coordinates (x A ) = (9, <j>) in such a way 
that n becomes a rotation field on the two-sphere. This 
means (cf. [16] or Appendix D) that there exists a coor- 
dinate system in which the following holds: 

ft = -Q, k e klm y l d m . (4.3) 

Here Vl k are components of a three-dimensional vector 
called angular velocity of the black hole, and y k are 
functions on S 2 created by restricting Cartesian coor- 
dinates on R 3 to a unit two-sphere. We can also set 
^-coordinate axis parallelly to angular velocity vector 
field. After a suitable rotation we have: (fi fc ) = (0, 0, 51), 
z = y 3 = cos 9, and: 

n = -a±. (4.4) 

Inserting this into (4.2) we obtain 

— -— f n A SW A = (4.5) 
87r Jov- 

where 

J=J X ~±- [ W v (4.6) 

is the z-component of the black hole angular momentum. 

Up to now we have used only the symmetry of con- 
formal structure carried by Hab- The symmetry of the 
metric itself implies that the conformal factor / is con- 
stant along the field n. This follows from the observation 
that the trace of the Killing equation implies vanishing 
of divergence of the field ft: 

= d A W&cth CD n A ) = n A \jdeth CD d A f, (4.7) 

where the fact that n is the symmetry field of the metric 
h has been used. Formula (4.4) implies that d v f = and 
the conformal factor / must be a function of the variable 
9 only. 
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It turns out that also its canonical conjugate k may be 
gauged in such a way that it is constant along the field 
n (see Appendix E for a proof). 3 

This result was obtained locally, or rather quasi-\oca\\y 
- i.e., from the analysis of the field on the horizon itself. 
However, the global theorems on the existence of station- 
ary solutions possessing a horizon, imply the so called 
zcroth law of thermodynamics of black holes (see [17]), 
according to which the surface gravity k must be constant 
along the horizon. But A := j s2 A is the area of the hori- 
zon S. Taking this into account and using (4.5), we derive 
from (4.2) the "first law of black holes thermodynamics": 

-sSM = -^-kSA + flSJ. (4.8) 

Contrary to the theory proposed by Wald and Iyer in [3] , 
the first law (4.8) is, in our approach, a simple conse- 
quence of the complete Hamiltonian formula (4.1), re- 
stricted to the stationary case. As illustrated by an 
example of the string dynamics, where formula (1.11) 
for virtual work was a consequence of the Hamiltonian 
formula (1.6), a similar "thermodynamics of boundary 
data" may be expected in any Hamiltonian field theory 
(see e.g., [2] for the corresponding analysis of the Maxwell 
electrodynamics) . Also a "Penrose-like" inequality (anal- 
ogous to (1.12) in the string theory) is satisfied as soon 
as the Hamiltonian is convex. We hope very much that 
the gravitational Penrose inequality can be proved along 
these lines. Preliminary results in this direction, based on 

3 In case £1 = 0, quantities k and / are arbitrary functions on S 2 . 



the analysis of the field Hamiltonian in linearized gravity 
(see [7]), are promising. 
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APPENDIX A: NOTATION AND PRELIMINARY 
FORMULAE 

A considerable simplification of the proofs is obtained 
if we use in a neighbourhood of S = {x 3 = const} a spe- 
cial coordinate system introduced in [14], which reduces 
the metric to the following special form: 





u a ua 




sM + m A nA 




n A 


9AB 


m A 




_ sM + m A n A 


m A 


(#) 


2 , A 

+ mm A 



(Al) 

Consequently, we have 









n A 

~w 


m A 
6 M 


s 

M 


<r = 


n A 


A M 


~ AB n A n B 
9 AT2 


i n n A m B +m A n B 
M 


* M 




s 

M 




-s^ 






(A2) 



where M > 0, s := sgn g 03 = ±1, g A s is the induced 
two- metric on surfaces {x° = const, x 3 = const} and 

~AB , ~AB 

g is its inverse (contravariant) metric. Both g and 
gAB are used to rise and lower indices A,B~ 1, 2 of the 
two- vectors n A and m A . 

We denote by A the two-dimensional volume form on 
each two-surface {x° = const, x 3 = const}: 

A := v/dctyAB. (A3) 
For any degeneracy field K of g a b the following object: 

A 

VK '-kW) 



is a scalar density on S. The vector density 

A = v K K = \(d Q -n A d A ), (A4) 

is well defined (i.e., coordinate- independent) and, obvi- 
ously, does not depend upon any choice of the field K. 
Hence, it is an intrinsic property of S. 

The external geometry of S is described in terms of 
the following tensor density: 

Q\(K) := -s {v K (V b K a - 5ZX7 C K C ) + 5 a b d c A c } (A5) 

which is fully analyzed in our previous paper [10]. 
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In our calculations wc shall use also quantities which 
are not geometric objects (are coordinate dependent). 
All of them drop out in the final result, where only well 
defined geometric objects remain. More precisely, we 
consider the following combination of the connection co- 
efficients: 

x (Tip - S 3 a T x px ) , (A6) 

and the two-dimensional inverse metric g AB rewritten in 
a "three-dimensional notation", where we put g := 0. 
The resulting matrix g does not define any tensor on S 
and satisfies the obvious identity: 

~ ac X a T^a rO 

g g C b = o b - k o b ■ 

Hence, the (0, 1, 2)-block of the contravariant metric (A2) 
may be rewritten as follows: 



—K a K b ~-(m a K» + m b K a ), 



(A7) 



where m a := g aE mj, so that m° := 0, and 



3 M = K». 
y M 

Using the above definition we may write that 

sQ\ =A (~g ac l cb - ±5 a b ^j + A a w b - S\A c w c 

+ A a Xb - S\A c X c , (A8) 
where l a b,w a are defined by (2.3) and (2.4) correspond- 
ingly, Xc '■— 2^ c ^ n ^ rom we a ^ so nave tnc 
following 

sQ a b = XS a b W c K c - XV b K a - S a b d c A c 
= -X5\{w c K c + l) + \(w b K a + f c l cb ) + S a b Xl 
= X~g ac l cb + A a w b - S\A c w c . (A9) 

Equation (A9) expresses the (2 + l)-decomposition of the 
three-dimensional density Q a b . As a consequence of (A8) 
and (A9) we obtain the following identity: 

sQ\ = sQ\ - l -\l5\ + A a Xb - S a b A c X c (A10) 

The detailed proof of these formulae is contained in paper 
[10], where we derive also the following equality: 

sQ af} g afi ,a = \(g be g c % d - \lg hc )9bc, a 
+ (A\f d + A c g bd - A d g cb )Al d9bc , a 
+ 2sQ 3 3 (d a \nM+jLm B n B a y (All) 



Replacing the partial derivative d a by the variation op- 
erator 8, we get an analogous formula 

sQ af} 5g aP = X(g be g cd l ed - \lg bc )8g bc 



+ (A b g cd + A c g bd - A d g cb )A 3 3d 8g bc 

1 , s 
M SM+ M- 



+ 2sQ 3 3 ( 4jSM + 47m B Sn B 



which may be further simplified to 



1 

<p V =sQ" , "dg ab - 

ia%bd I \b~ad \d~ w ab\ 



sQriQw =sQ ab 6g ab ~ ^-Xl~g ab 8g ab - XIS log M 



+ X d (A a g bd + A b ~g ad - A d ~g ab ) 6g ab . (A12) 
Moreover, we need the following identities from [10]: 

crv = - sX ( l + 2Rd ( w d + x^) , (Ai3) 

and 

Xl = -X + d A (n A X). (A14) 



APPENDIX B: PROOFS OF THE GENERATING 
FORMULAE FOR DYNAMICS OF 
GRAVITATIONAL FIELD 

Dynamics of gravitational field is derived from the 
principle of the least action 8A = 0, where the action 
of gravitational field is defined as the integral of Hilbcrt 
Lagrangian: 



1 

T6tt" 



R. 



(Bl) 



The method proposed by one of us in papers [1, 2] leads 
to Einstein equations written in the following form: 



SL = d K (n^SA;„) 



(B2) 



where 



\9\ 9 



fJbV 



1 

16tt" 

v x — fi x r K 



As soon as we choose a (3+l)-decomposition of the space- 
time M , our field theory will be converted into a Hamil- 
tonian system, with the space of Cauchy data on each of 
the three-dimensional surfaces playing role of an infinite- 
dimensional phase space. Let us choose coordinate sys- 
tem adapted to this (3+l)-decomposition. This means 
that the time variable t = x° is constant on three- 
dimensional surfaces of this foliation. We assume that 
these surfaces are spacelike. To obtain Hamiltonian for- 
mulation of our theory we shall simply integrate equation 
(B2) over such a Cauchy surface C t C M and then per- 
form Legendre transformation between time derivatives 
and corresponding momenta. 
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We consider the case of an asymptotically flat space- 
time and assume that also leaves C t of our (3+1)- 
decomposition are asymptotically flat at infinity. To keep 
control over two-dimensional surface integrals at spatial 
infinity, we first describe dynamics of our "matter + grav- 
ity" system within a finite volume V. Integration of (B2) 
over V yields: 



5 [ L = [ d K (tt^SA^) = 
Jv Jv 



n^SA^, (B3) 



where "dot" denotes time derivative (the two- 
dimensional divergence 8b (tt^SA^ vanishes when 
integrated over dV). 

In [1] the Legendre transformation between time 
derivatives and corresponding momenta was performed 
in case of a non-degenerate (one-timelike, two-spacelike) 
surface S. Here we generalize this method to the case of 
a wave front. The first step in this construction consists 
in observation that, due to metricity of the connection T, 
the following identity holds: 



^6A% = -^-g kl 6P kl +d k (n 00 5 r " 



, (B4) 



where P kl denotes the external curvature of S written in 
the ADM form. A simple proof of this formula is also 
contained in paper [10]. 

On the other hand, direct calculations of the variation 
of the quantity Q^ v given by (A6) lead to the following 
reduction of the boundary term tt^SA^ = ir^ v 8A 3 : 



1 



^5Al„ = -—g^5Q^. 



167T' 



(B5) 



Skipping the two-dimensional divergencies which van- 
ish after integration and using (B4) and (B5), we may 
rewrite the right-hand side of (B3) in the following way: 

/ (n^6A%y+ [ n^8A^ 

JV J dV 

--ib/ v <*^>" + jU*"'(£))' 

g^SQ^- (B6) 



16tt 



Now, we perform the Legendre transformation both in 
the volume: 

(9kiSP kl )' = (g k i5P kl - P kl 8g kl ) + S (g kl P kl ) (B7) 
and on the boundary: 

^03 



(B8) 



.03 



In paper [10] the following formula has been proved: 



{to**) 



16-7T / 7T 
dV 



.00 



T 03 



.00 



= {d k (V\9\ .9°^)) + 2VW\R°o} ■ (B9) 

Hence, generating formula (B3) takes the form: 



1 

167T Jgy 

7T M S [ — ) - ( — ; Sir' 



03 \ / n 03- 



dV 
1 



7f 



00 



.00 



+ 



(BIO) 



Using the form of the metric (Al) and (A2) we express 
tt^ in terms of the metric. Denoting 



a := \ogN 



(Bll) 



we obtain 



^r00 X 

IT O 

<>v \ V " 
s 



7i~ 



00 



.00 



Sir 



00 



16tt 



<)V 



2XSa 



- 2aSX) 



(B12) 



- — f (kSlogM - (log My 8 X 

lOTT JdV V 

Similarly, we rearrange the boundary term: 
g^QT = S (g^Gr) - Q> lv 8g IJ , v . 

For this purpose we use the formula (A12) and (A13). 
Therefore, our generating formula takes the following 
form: 

^j v MR-^) = -^j v Q\ 

(k8a - aSX") + JL J Q ab (K)8g ab 



iv 
-f 

87r JdV 

s 

16tt ./,,, 



IdV 

l^9 ab 8gab - AM log M - XSlogM 



+(log My 8 X + Xd (A a g bd + A b g ad - A d g ab ) Sg ab 

it)7r JdV 

+ ^ J gy VW\ (g 3 ^l g°^l) ■ (Bi3) 
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Now, we simplify the following boundary term (in square 
brackets in the above formula): 

-^ ab Sg ab - XIS log M - \6\ogM + (log My SX 

(B14) 



+ Xd (A a ~g bd + A b ~g ad -A d ~g ab )8 9ab . 



Using identity (A14) (see e.g., [14]), A b 8g bc = -g bc 5A b , 
and skipping two-dimensional divergencies, we obtain 
that the above expression takes the following form: 

X(d A \ogM)8n A - d A (n A X)S log M 

+ n A d A (\ogM)SX + (d A n A )SX - X(d A log X)Sn A . 

(B15) 

Further simplifications can be made after integration over 
dV: 

f (d A (\ogM)SA A - (d A A A )S\ogM) 

JdV 

= -6 ( ((d A A A ) log M) , (B16) 
Jdv 

[ {(d A n A )SX - (d A X)Sn A ) 
Jdv 

= f (d A n A )5X + XS(d A n A )) 
Jdv 

= 5 [ X(d A n A ). (B17) 

JdV 



Finally, expression (B14) takes the following form (mod- 
ulo two-dimensional divergencies): 

5 ((d A A A ) log M + X8 A n A ) , 

and the formula (B13) reads as 



+ ^-J (^a - a5X) + -L j Q"<'( K )6 !hlh 



8tt 



16tt 



dV 



+ 



+ 



^-5 j ((d A A A )\ogM + Xd A n A ) 

1D7T Jgy 

f X (2K d (w d + X d) + I) 

51" JdV 

+ ii a / ) ,^(" 3 " r8 --»°" r w- 

Now, we simplify the expression: 



(B18) 



S I {s(d A A A ) logM + sXd A n A 
hv 

+ sX {2K d (w d + Xd ) + 
+ VW\(g 3 ^l-g^Tl)}. 



Using 



TL - —m% 

3a M 



M 



-M, 



(B19) 



(B20) 



and taking into account that K^T = —w a (2.4), 

Kp^ab = l ab (^-3) we can pass to the metric derivatives 
contained in T's. This way (B19) reduces (modulo two- 
dimensional divergencies) to 

-2s [ lSX + 2s5 [ X(l-n A w A ). (B21) 

JdV JdV 

Therefore, we obtain the following boundary formula: 
+ 77- I (X5a-a5x) 

+ \\hj dv ( Qab ( K ^ b - 2slsx ) 

+ ^-S [ X(l - n A w A ). (B22) 
87r Jdv 

Using identities Q ab (K)5g ab = 2Q° A (K)Sn A - 
Q AB (K)5g AB and sXn A w A = Q° A (K)n A we get: 

+ — / (x5a-a5x) 

8?T Jgy V ) 

[ (Qab (K)5~g AB + 2n A SQ° A (K) - 2sX6l) . 

ibn JdV 

(B23) 

Using equation (A9) we express Q in terms of indepen- 
dent objects l a b and w a , and finally obtain: 



- sn = J (P kl Sg kl - g kl SP kl ) 



+ — / (X6a - aSX) 
8tt Jav 



+ tJ^/ (U AB SgAB-2(w 5A°-A A Sw A )) 7 

(B24) 



where 



1 



»7T Jy 87T Jgy 87T Jgy 



APPENDIX C: DERIVATION OF GENERATING 
FORMULA FROM THE NON-DEGENERATE 
CASE 

Consider a one-parameter family of hypersurfaces 
S e :— {r e := r — set — const} parameterized by a real e, 
such that for e = we have So = S and for each e ^ 
the induced three-metric on S e is non-degenerate. 

Take the external curvature tensor Q ab of S e and the 
three-dimensional contravariant metric g ab , inverse to 



10 



g ab . Define 



Q ■= , 

Qa := Q°A, 
i 

~>AB 



Q := Qcd~9^~9 



CA-DB 



(CI) 
(C2) 
(C3) 



Observe that these are two-dimensional objects defined 
on dV: Q is a scalar density, Q A is a covector density 
± 

and Q AB is a symmetric tensor density. The following 
identity (a homogeneous generating formula for the field 
dynamics) was proved in [2] for a timelike boundary and 
generalized in [18] for any non-degenerate (e.g., a space- 
like) boundary: 

= J (P kl 5g kl - g kl SP kl ) + 2 J (k6a - a<5A) 

+ / (Q AB S 9 ab - 2n A SQ A + 2nSQ), (C4) 

JdV 

where 



n := 



„A %AB n 

n -—9 9ob- 



The boundary objects a, Q, Qa and Q AB diverge under 
the limit S e — > S - or, equivalcntly, g 33 — > 0. We are 
going to rearrange them in such a way that the resulting 
objects do behave well under this limit. 

For this purpose three-dimensional metric component 
g 00 defining the lapse function n, can be expressed by the 
following components of the spacetime metric: 



~oo = g oo. (.9° 3 ) 2 _9 m 9 33 -(g° 3 ) 2 



.,33 



~33 



(C5) 



Four-dimensional component g 33 of the inverse metric 
can be expressed in terms of three-dimensional metric 
component g 33 inverse to metric gu induced on V: 



9 33 = ~9 33 + 



(g 03 ) 2 



(C6) 



Then 



00-33 

~oo _ 9 9 

* „33 ' 



(C7) 



Restricting formula (Al) to the hypersurface {x° 
const} we obtain: 



9kl 



9AB 



m A 



m A 



(f) 2 + 



m A mA 



Consequently, its three-dimensional inverse equals: 



((§) 2 +m A m A )~g AB 


A 


-m A 


1 



In this notation g 00 = —1/N 2 , hence we obtain that 
g 00 = -l/(M 2 g 33 ) and the lapse n equals 



n = M^/g 33 . 



(C8) 



Using definition of the hyperbolic angle a — arsinh(q) , 

30 

where q = — ,. __ we obtain: 



33 



6a 



6q 



Taking into account that 



5(logg). (C9) 



log q = log 5 03 - - (log |.g 00 | + logg 33 ), 
we obtain: 

q 6a = 6 log N - 6 log (M^fg* 3 ) 
Vl + 5 2 V ' 

= 5a — <51ogn, 

where 

o:= -logVlP°l =logN 



(CIO) 



(Cll) 
(C12) 



is a regular part of a. The second term Slogn diverges 
in the limit, but cancels the divergent part of Q, what 
will be proved in the sequel. 
Moreover, we have: 



XI, 



(C13) 



nQ = n 2 Q 00 = XMT 3 AB ~g AB 
hence the following formula holds: 

2n5Q = 2nQ6 {log Q) = 2X16 (log XI - logn) . (C14) 
Continuing, we have: 

Qa = Q° A = g ob Q b A 

= AM (> - ^f) (r 3 bA g bA T 3 cd ~f d ) . (CIS) 

Using the following formula: 

-i(io gff 33 ), a = r 3 3 + -Lr 3 b .g 3h , (Ci6) 

z 9 



we obtain 



Q°a = A (r 3 A - m B Y 3 BA ) + -X (log 



„33\ 



■A ' 



(C17) 
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Now, we define the quantity: 



(C18) 



Taking into account that 

TL - m B T% a ^W a + (log M), 0) (C19) 
it follows that Q° A may be written in the following form: 

Q% = Arc a + A(logM), A + ^(log ff 33 ), A . (C20) 



hence the contraction Q AB g AB gives us 



Q AB 9AB = -XI - 2Atr 

+ 2\n A K> A - 2X(d - n A d A )logM 
-X(d -n A d A )(\ogg 33 ), 



and the traceless part of Q AB takes the form: 



Similarly, Q AB takes the following form: 

± 

Q AB = Qcd~9 ca § db 

= \M~f A ~g DB (T* CD - gcD Tl f Ff) , 



(C21) 



Q AB - \ Q CD 9cd9 AB = A (l AB - \fg AB 



The above results may be gathered as follows: 



(C22) 



2\8a = v/l + g2 (2\5a - 2X6 log nj 



-2aSX = 



q 



(-2d aSX + 2d logn5X) , 



2nSQ =25 (XI) - 2AW(logn), 
~2n A 5Q A = - 2n A 5(X\v A ) - 2n A S(X(logn) A ), 

- [I + 2tv - 2n A rv A + 2(d - n A d A )(\ogn)] SX, 



Q AB g AB S\ogX 



1 J- 



Q AB - ^ Q CD 9cd~9 AB ) 59ab =X ( l AB - ^lg AB ) 5g AB . 



<AB 1 i~AB 



Finally, we express quantities appearing in boundary for- 
mula (C4) in terms of the above objects. Taking into 
account identity (A14) and omitting two-dimensional di- 
vergencies we obtain 



= J (P kl 5g kl - g kl 5P kl ) 



(C23a) 



+ 2 I y/l + q \ x5a - a8X) + 25 f XI 
hv <7 JdV 



(C23b) 



+ I (Xl AB 5g AB -2Xn A 5ro A -2w 5X) 

JdV 

(C23c) 

+ 2 [ \X51ogn- ^ 1 + q2 XSlogn ) (C23d) 
Jdv \ <1 I 



dv \ q 



-(\ogn)5X — (log n) 5 X 



(C23e) 



In the null limit e — > 0, we have that a = a and ro a = w a . 
Moreover, q — > oo and, whence, y/l + q 2 /q — > 1. We 
shall prove that the last two boundary terms (C23d) and 
(C23c) vanish in this limit. Indeed, we have 



lim 



I " 1 I Slogn= — <V|ff 33 | 



(C24) 



A similar formula with 5 replaced by do is also true. But 
Sg 33 -> and g 33 -> on dV. This implies: 



lim 

£^0 



<51ogn 



=0, 



lim 



1 <9 (logn) 



=0. 



9V 
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Hence, boundary formula (C23) takes the following form: 
0=J v {P kl Sg kl -g kl SP kl ) 

+ 2 1 (XSa - aSX) + 25 j XI 

JdV JdV 

+ [ (\l AB 5g AB -2\n A 5w A -2w 5\). (C25) 

JdV 

Moving 26 J gv XI to the left-hand side, we obtain formula 
(2.1). 



APPENDIX D: DERIVATION OF THE 
FORMULA (4.4) 

Being a difference between two symmetry fields: the 
four-dimensional Killing field do and the null field K on 
S, the field n — (n A ) is a symmetry field of the two- 
metric gAB and, whence, satisfies the two-dimensional 
Killing equation: 



n>A\\B + n B \\A = 0. 



(Dl) 



As was already discussed, we may choose a coordinate 
system such that gAB = fh>AB, and Hab is a standard 
two-sphere metric. The field n is also the symmetry field 
of conformal structure given by the metric Iiab-, hence it 
has to fulfill the equation: 



n A]\B+ n B\\A- h ABn C r lc = 0, 



(D2) 



where || denotes the two-dimensional derivative with re- 
spect to metric h on the two-sphere. Therefore, n belongs 
to the six-dimensional space of conformal fields on two- 
sphere and it is of the following form: 



„A _ AB , f. -uAB 

n — s v g + v b n , 



(D3) 



12 1 

where v, v arc dipolc functions, i.e., of the form: v= aik l 7 

2 

v= b; t k l \ where k l are coordinates of the unit vector: 

k 1 = sin 9 cos ip, 
k 2 = sin 9 sin ip, 
k 3 = cos 9. 

Taking all this into account we write the field n in the 
following form: 



n 



We will show that there exists a coordinate system in 
which n may be written as 

n A d A = 



Proof. Suppose that for 9 = component n vanishes. 
Hence a\ = b 2 i a 2 = — b\. Then the components of n 
read as 

n e = (cos 9 — l)(ai sin tp + b\ cos ip) — b 3 sinO, 
■rif = (cos 9 — l)(bi sin p — a\ cos tp) + a 3 . 

We can rotate the coordinate system in such a way that 
the components n and will take the form: 

n e = (cos 9 — I) a sin p — b sin 9, 

n v = (cos9 — 1)(— acosp) / - + c, 

sine* 

a, b, c are some new parameters. Because S is a horizon, 
we have: 

(\n A ). A = 0, 

and in our coordinate system A = / sin 6, hence we have 
the following equation: 

= {/sin 6* ((cos6> — l)asin<£ — 6sin6 1 )} e 

+ |/sin6»^(cos^-l)(-acos<^)^^+c^| . (D4) 

Integrating the above equation over ip and omitting van- 
ishing integral J (\n ip ) tip dtp we obtain 

d f 27T 

— / fsin9 ((cos9 — l)asmtp — bsinff) dip = 0, 



hence 

cos 6»-l r27r 



(D5) 



sin 9 



27T 



/ sin ipdip + b fchp = 0. (D6) 



Because the first summand vanishes for 9 — > 0, we have: 

r-2-K 

2nb f(9 = 0)=b fd<p = 0. 
Jo 

This straightforwardly implies 6 = 0, because the con- 
formal factor / is positive. Hence, n is of the following 
form: 

n e = (cos 9 — l)asimp, 

n v = (cos 9 — 1)(— a cosy) + c. 

sm9 

Let us write n in stereographic coordinates on the plane 
(x, y) intersecting our two-sphere along equator, where 
coordinates (x, y) read as: 



dip 



x = ctg^ cos ft 
V = ctg-smy. 



(D7) 
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Then: 



d_ 
~d~9 



1 cos ip d 1 sintp d 



2 sin 2 fete 2 sin 2 § cty ' 



9 6> . 9 9 d 

— — = — ctg- simp- — h ctg- cos <p-^~, 



<9</5 



9a; 



dy 



and n A &A is of the following form: 



4 ^ <9 
n Oa = a— + c 
dy 



— +x— 

dx dy 



(D8) 



Let us consider two cases: c = and c ^ 0: 



• c = 

In this case equation (D4) takes the following form: 



(/sin^cosfl — l)asin ip) g 

+ I / sin 9(cos 9 — 1) (—a cos ip) 



sm( 



0. (D9) 



In stereographic coordinates (D7) the above equa- 
tion reads: 

Oogf) y = z - 2 ^ 2 , (D10) 
' y 1 + x A + y z 

and implies 

f(x,y) = C(x)(l+x 2 +y 2 ), (Dll) 

where C(x) is an arbitrary function of x. Keep x 
constant and pass to the limit y — > oo. If C(x) 7^ 
then / — > 00, otherwise / = 0. But the confor- 
mal factor / must be finite and different from zero. 
Hence, the case c = is incompatible with the 
properties of /. 

• 

In that case we may write 



n A d A 



hi) 



a\ d 
c) dy 



dx 



Putting x = x + -, we obtain: 



d 



n A d A = c [ x^- - y-^z 
dy dx 



(D12) 



(D13) 



hence n A dA is a field of rotations and can be writ- 
ten as follows: 



n A d A = -H 



d_ 

dip 



(D14) 



A different proof of this fact is given in the paper [16]. 



APPENDIX E: SURFACE GRAVITY k IS 
CONSTANT ALONG ft 

We are going to perform a gauge transformation of k, 
such that the resulting quantity k remains constant along 
n. For this purpose we leave the spacelike coordinates 
unchanged (x k = x k ), whereas the time coordinate is 
translated by a constant which depends upon them. This 
implies the following transformation on the surface S, 
where x 3 is constant: 

x°=x° + a(x A ). (El) 

Then the field K transforms as K = cK, where 

c= (1 , (E2) 

and quantities w a transform as follows: 

w a = w a + d a (log c) (E3) 

(transformation law for objects w a is given in [14]). Af- 
ter this transformation the surface gravity k takes the 
following value: 

k= -K a w a = -cK a (w a + ^j =ck-K{c). (E4) 

Using K = do — n A dA and formula (D14) for n A we 
obtain 



n 9c . n 

il— CK + K = 0. 

dtp 



Consequently, 



Denoting 



we obtain 



F(9,ip) = J e(" n 'f^dip 



C - 1 = -?^(logF) 
k dip 



(E5) 



(E6) 



(E7) 



(E8) 



We can compare this with the form of c 1 implied by 
(E2) 



! = 1 



which leads to the following equation: 

da 1 IS, r 
— = 1 log F. 

dip k dip 



rj Solutions of the above equation reads as 



ip 1 
a= £ + ^logF, 



(E9) 



(E10) 



(Ell) 
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where 



F{<p) = Ci QT e ( " a 1 /o" K ( s ) ds W + C*2^) , 

and Ci, C2 are integration constants. Applying period- 
icity conditions: 

a(0) = a(27r), 

A (logF)(0 ) = ^(i ogi ,)( 27r) 



we obtain 



- n . f(o) 

AC — — — log 



2tt b F(27r)' 



(E12) 



where the values of F(0) and F(2n) remain to be deter- 
mined. Denote by f(tp) the expression: 

/(¥>):= ^< s)ds) du. (E13) 

Jo 

Therefore F and log F are of the following form: 

F = C 1 (f + C 2 ), 
logF-logC 1 +log(/ + C 2 ). 



The above equations imply 



/'(0) 



/(0) + C 3 /(27r) + C 3 



and /(0) = 1, (E14) 



hence 



C 3 



/(27T) 
/'(27T)-!' 



and we have 



- n . c 3 

K — — log 



(E15) 



(E16) 



2tt fa /(2^) + C7 3 ' 

Finally, we obtain 

f2 1 f 27T 

k = -— log/'(27r) = — / /c(s)ds ee const. (E17) 

27T 27T J 

Hence, we have performed such a gauge transformations 
from k to k, that the new quantity k is constant along 
parallels of the sphere S 2 . 
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